Interaction-induced correlations and non-Markovianity of quantum dynamics 
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We investigate the conditions under which the trace distance between two different states of a 
given open system increases in time due to the interaction with an environment, therefore signalling 
non-Markovianity. We find that the finite-time difference in trace distance is bounded by two sharply 
defined quantities that are strictly linked to the occurrence of system-environment correlations cre- 
ated throughout their interaction and affecting the subsequent evolution of the system. This allows 
to shed light on the origin of non-Markovian behaviours in quantum dynamics. We best illustrate our 
findings by tackling two physically relevant examples: a non-Markovian dephasing mechanism that 
has been the focus of a recent experimental endeavour and the open-system dynamics experienced 
by a spin connected to a finite-size quantum spin chain. 
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Recently great attention has been paid to the develop- 
ment of a more general understanding of the dynamics of 
open quantum systems, in order to deal with the occur- 
rence of memory effects [TlU7|. In particular, different 
definitions of quantum non-Markovianity have been the- 
oretically introduced [18-21 and, in some cases, experi- 
mentally investigated [22 J 2^. Nevertheless, the physical 
reasons ruling whether an open quantum system exhibits 
a Markovian or a non-Markovian dynamics have still to 
be fully clarified. 

A widely accepted view is that, in Markovian dynam- 
ics, the correlations between the open system and its en- 
vironment as well as the changes in the environmental 
state due to the interaction do not have a significant in- 
fluence on the subsequent evolution of the open system. 
This picture is often introduced relying on qualitative 
considerations, possibly assuming that the total state at 
time t can be effectively represented as a product state 
between the state of the open system at the time t and 
a fixed state of the environment [TJ . It is worth not- 
ing how the same assumption also lies at the foundations 
of the quantum regression hypothesis 27-30 . System- 
environment correlations induced by the interaction and 
changes in the state of the environment are thus thought 
to be at the basis of non-Markovian dynamics. 

In this paper, we show how this relationship can be 
formulated in a quantitative way using the properties 
of the trace distance [31] . whose time evolution can be 
used to characterise the dynamics of an open quantum 
system that evolves starting from two different initial 
states 13U [35] . Any change of the trace distance be- 
tween states of an open system can be interpreted as an 
exchange of information with the environment that af- 
fects it: a non-monotonic behaviour of the trace distance 
witnesses the fact that some information previously lost 
by the open system can affect it back again, thus induc- 
ing memory effects in its evolution. In view of this in- 
terpretation, non-Markovian dynamics can be identified 
with those dynamics that show an increase of the trace 



distance at some intervals of time [HI [33] . 



Here, starting from the analysis reported in Ref. |15j . 
we introduce an upper and a lower bound to the vari- 
ation of the trace distance at finite time intervals. 
The bounds express quantitatively the influence of the 
system-environment correlations and the changes in the 
state of the environment at a time t on the subsequent 
dynamics of the open system. They thus allow to es- 
timate how system-environment correlations, as well as 
evolution of the environment, account for the Markovian 
or non-Markovian nature of open system's dynamics. In 
particular our lower bound provides a sufficient condi- 
tion for the onset of non-Markovianity. We apply our 
analysis to two physical examples: the experimental set- 
ting considered in Ref. [22] that describes the transition 
from Markovian to non-Markovian dephasing on a qubit 
and the energy-non-conserving open-system dynamics of 
a spin-1/2 particles that is coupled to a finite-size quan- 
tum spin chain [9] I34j . 



The remainder of this paper is organized as follows. 
In Sec. [I] we first introduce an upper and a lower bound 
to the variation of the trace distance, which hold under 
very general conditions. We further show how, as a conse- 
quence, the strength of the effects of system-environment 
correlations and environmental evolution can determine 
the non-Markovian character of a given dynamics. In 
Sec. [II] we apply our general analysis to two physically 
relevant examples. First, we address the single-qubit 
pure-dephasing mechanism exploited in Ref. [22) to inves- 
tigate experimentally the transition between Markovian 
and non-Markovian dynamics. Second, we study the case 
of a single spin interacting with a finite-size spin environ- 
ment embodied by a quantum spin chain [9, 34! . Finally 
Sec. Mil is devoted to conclusions and final remarks. 
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I. ROLE OF SYSTEM-ENVIRONMENT 
CORRELATIONS AND ENVIRONMENTAL 
EVOLUTION IN THE DYNAMICS OF OPEN 
QUANTUM SYSTEMS 

A. Upper and lower bound to the increase of the 
trace distance on a finite time interval 

We now highlight the relevance of the correlations be- 
tween system and environment, as well as of the changes 
in the state of the environment, in determining the vari- 
ation of the trace distance among system states. As dis- 
cussed above, knowledge of the time dependence of this 
quantity allows to assess the Markovianity properties of 
the time evolution, and is therefore of primary impor- 
tance. To this aim let us start from the decomposition 
of any joint system-environment state pse as [35] 



Pse = PS® PE + XSE, 



(1) 



where ps (pe) is the reduced state of the system S 
(the environment E) and xse (such that tr^ [xse] = 
trs [xse] = 0) accounts for the total correlations between 
the open system and the environment, as also shown by 
the relation ||xse|| = 2D(ps <S> Pe,Pse)- The total S-E 
system is usually assumed to be closed, so that its evo- 
lution is provided by a one parameter group of unitary 



operators {Ut} t>0 , with to = initial time. Given the 
initial total state pse(0), the total state at a time t is 
PSE{t) = U t psE(0)Uj and, as a consequence, the state 
at time t + 1' can be inferred from the state at the time 
t through the relation 



PSE(t + t') = U V ,tPSE{t)Ul. 



(2) 



where U t >,t = U t+ t>Uj. 

Using Eq. ([I]) twice, we can simply express the differ- 
ence in the total states at time t originating from different 
initial conditions as follows: 



P SE (t) ~ P SE (t) = (Ps(t) 

+Ps(t)®(PE(t)-pUt)) 



Ps(t))®P E (t) (3) 

(XSE(t)~X 2 SE(t))- 



An equivalent relation is obtained by exchanging the role 
of labels 1 and 2. The difference between the total states 
can thus be split in two contributions, one depending on 
the difference between the states of the reduced systems 
and the other, made up of the last two contributions at 
the right hand side (r.h.s.) of Eq. (jsj), given by the com- 
parison between the reduced environmental states and 
between the correlations. We thus identify the two quan- 
tities 



and 



F^,t,p x ' 2 )=D{tv E [Ut>M{t)®P E {t))Ui.l\ Me [Ut>APs(t)®PE(t))ui,t]) (4) 

Uv,t (xsE(t)-X 2 sE(t))ul 



B(t',t,p^ 2 ) = -|| tr B \U v , t (p 2 s (t) <8 (p E (t) - P 2 E {t))) Uy 



tr, 



(5) 



where we have introduced the trace norm of an element 
a of the set T of linear trace class operators as ||cr|| = 
Tr[V a^a] and the trace distance between two statistical 
operators p 1,2 

D(p\p 2 ) = l\\p 1 -P 2 \\ = \j2\^l ( 6 ) 

fc 

with Qk the eigenvalues of the self-adjoint traceless oper- 
ator p 1 — p 2 and < D{p 1 ,p 2 ) < 1. D(p 1 ,p 2 ) quantifies 
the distinguishability [31] between two given states. If a 
system is prepared in one of the two states p 1 or p 2 with 
probability 1/2 each, an observer can design an optimal 
strategy to guess the preparation with success probability 
given by [1 + Dip 1 , p 2 )}/2. 

The quantity F(i',i,p 1,2 ) describes how the distin- 
guishability between the reduced states would evolve at a 
time t+tf if the two total states at the time t were product 
states, with the same environmental state, thus building 
on the first term at the r.h.s. of Eq. ([3]). In a complemen- 
tary way, the quantity B(t', t, p 1,2 ) keeps track of the ef- 



fects of correlations and differences in the environmental 
states at a time t on the subsequent dynamics of the open 
system. Thanks to the contractivity of the trace norm 
of a self-adjoint operator under the action of any com- 
pletely positive trace-preserving (CPT) linear map |36) . 
it is straightforward to show that B(t' ,t, p 1 ' 2 ) G [0,2]. 
Moreover, it can be null even if the total states at time 
t are not product states, so that the correlations, despite 
being present, do not have any influence on the evolution 
of the trace distance between two reduced states of S. 

In what follows, we will study the evolution of the dis- 
tinguishability between couples of reduced states, thus 
describing the information flow between the open sys- 
tem and its environment [35J [33] . The trace distance 
between two reduced states at time t will be indicated as 

D^p^^DipUt)^^)), (7) 

where Pg(t) = tr B [[/ t p^ B ((])[//] , j = 1,2. Our analysis 
focuses on the variation of the trace distance at finite 
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time intervals 

AD(f , i, p 1 ' 2 ) = D(t + If, p 1 ' 2 ) - D(t, p 1 ' 2 ). (8) 

Using Eq. ^ and ([3]), the trace distance between two 
reduced states at time t + t' can be expressed as the sum 
of two different contributions, which reflect the decom- 
position at the r.h.s. of Eq. p|. Then, as the trace norm 



is unitarily invariant and the partial trace is a CPT map, 
the contractivity of the trace norm under CPT maps as 
well as the triangular inequality for trace norm 

IH-lk'HI < < H + ||<7'|| Wa,cr'eT, (9) 

directly lead to the following bounds to the variation of 
the trace distance on finite time intervals 



B(t', t, p 1 - 2 ) ~ F(t', t, p 12 ) - D(t, p 1 ' 2 ) < AD(t', t, p 1 - 2 ) < B(t', t, p 1 ' 2 ) + F(t', t, p 12 ) - D(t, p 1 ' 2 ). (10) 



Eq. ( 10 ) holds in complete generality, the only require- 



ment being that one takes into account the full unitary 
evolution Ut, which could well be time inhomogeneous. 
The upper bound in Eq. ( 10 ) shows that 



B(lf,t ) p 1 ' 2 )>D(t,p 1 > 3 )-F(l/ ) t,p 1 < 2 ) (11) 

is a necessary condition to the increase of the trace dis- 
tance within the time interval [t, t + t'] . It is important 
to note that the r.h.s. of Eq. (Ill is positive as a con- 



sequence of the contractivity of the trace norm. In fact, 
F(t',t,p 1,2 ) can be written as [15] 

F(t',t,p 1 ' 2 )=D{$ t ,, t [pUt)] ,*v,t [P 2 s(t)]) , (12) 
where 



*t',i [ps] = trjs U t ,.t{ Ps ® pUtWl 



(13) 



for any t and t' , so that $ t < t is a CPT map pQ and 
Fff.t.p 1 ' 2 ) < D&p 1 ' 2 ), with i^O^p 1 * 2 ) = Dit^p 1 ' 2 ). 
An increase of the distinguishability between reduced 
states needs the effects of system-environment correla- 
tions to prevail on the contraction of the trace distance 
due to the reduced CPT map <&t',t obtained from a total 
product state at time t. Notice that in the limit t' — > 
the upper bound in Eq. ( 10 ) leads to the bound found 
in Ref. |15) . Finally, by applying again the triangular 
inequality and using the contractivity of the trace norm 
under CPT maps on Eq. ( 10 1, one gets the weaker upper 
bound 

AD(t', t, p 1 - 2 ) < D{p% E {t), pKt) ® p l E (t)) (14) 
+D(p 2 SE (t),p 2 (t) ® p 2 E (t)) + D(p E (t),p 2 E (t)). 

This confirms that an increase of the trace distance in 
the time interval [t + t',t] calls for system-environment 
correlations in at least one of the two total states at time 
t, or for different environmental states p E (t) and p E (t). 
The inequality in Eq. ( 15 ) was first derived in Ref. |32j 
taking t = as the initial time of the dynamics and 
pointing out that an increase of the trace distance above 
its initial value witnesses initial correlations or different 
initial environmental states. Indeed, a non-monotonic 
temporal behavior of the trace distance can be read as 
an increase with respect to its value at a previous time. 



Up to now we have seen how system-environment cor- 
relations can be in general compatible also with a de- 
crease of the distinguishability between reduced states. 
But to what extent is this the case? Can the effects 
of correlations and environmental evolution be arbitrar- 
ily strong without inducing any trace-distance increase? 
This question is answered by virtue of the lower bound 
in Eq. (10). A sufficient condition to have an increase of 



the trace distance in the time interval [t, t + t'] is that 
the effects at the time t + t' of the system-environment 
correlations and environmental states at the time t are 
strong enough to satisfy 

B(t / ,t,p 1 < 2 )>D(t,p 1 > 2 ) + F(t / ,t,p 1 > 2 ). (15) 

Let us emphasize how the fulfillment of this condition 
implies the occurrence, with certainty, of an increase of 
the trace distance, which is quite a remarkable result. 



B. System-environment correlations and quantum 
non-Markovianity 

In virtue of the bounds to the variation of the trace 
distance discussed in the previous Subsection, we can 
now show that system-environment correlations and 
changes in the environmental states actually determine 
the Markovian/non-Markovian character of a dynamics 
according to the definition given in Ref. [15] . 

Let us assume a product initial state pse (0) = ps (0) <g> 
Pe(0) with a set environmental state p^(0), which implies 
the existence of a well-defined reduced dynamics on the 
whole set of statistical operators of the open system. In 
fact, the open system's dynamics can then be described 
via a family of CPT maps {A(i)} t>0 , with [T] 



A(t) [p s ] =tr E \u t (p s ® P E (0))U} 



(16) 



so that ps(t) — A(t)ps(0). Non-Markovianity of quan- 
tum dynamics has been defined and quantified in terms 
of the various properties of this family of CPT maps. In 
particular, the measure of non-Markovianity A/"(A) intro- 
duced in [TH] can be expressed as 

AT (A) =maxVA£>(6 fc -a fc ,a fc ,p 1 ' 2 ) ) (17) 

n l,2 ( * 



4 



where (au, 6fe) are the time intervals where the trace dis- 
tance D(t,p ' ) increases, and the maximum over all 
pairs of initial reduced states is taken. A/"(A) quanti- 
fies the total amount of information that flows to the 
open system, as witnessed by the trace distance, and 
the relevance of memory effects on the reduced dynam- 
ics. From Eq. (17 1, it is clear that a quantum dynamics 



is non-Markovian if and only if there is a time interval 
[t,t + t'] and a pair of initial states Pg 2 (0) of the sys- 
tem such that ADit^t^ 1 ' 2 ) > 0. As we have assumed 
an initial total product state, with fixed environmental 
state, the dependence on p 1 ' 2 should be intended from 
now on as referred to the two initial reduced states. In 
addition, possible system-environment correlations in the 
total states, as well as differences between the environ- 
mental states, at a time t have entirely to be ascribed to 
the system-environment interaction up to time t. 



Eqs. (11) and (15) provide a general reference scale 



that relates the Markovian or non-Markovian nature 
of a given dynamics to the relevance of the quantity 
Btf^p 1 ' 2 ) with respect to D^^p 1 ' 2 ) and F(t' ,t, p 1 ' 2 ). 
This relation is schematically depicted in Fig. [T] where 
the horizontal axis represents the possible values of 
B(t' ,t, p 1 ' 2 ) together with the two reference values given 
by D{t,p x ' 2 ) ± F(t' '^p 1 ' 2 ). The line above (below) the 
upper (lower) threshold denotes the region where non- 
Markovianity (Markovianity) is enforced. For the inter- 
mediate region of width 2F(t' ,t, p 1,2 ) no general state- 
ment can be made from the two bounds. Indeed, even 
if system-environment correlations and evolution of the 
environment due to the interaction affect the dynam- 
ics of the open system, this does not guarantee that 
the dynamics is non-Markovian. But the reduced dy- 
namics is surely non-Markovian if the effects of system- 
environment correlations and environmental evolution, 
as quantified by B{t' ,t, p 1 ' 2 ), are strong enough to ex- 
ceed the upper threshold in Fig. [T] 

In addition, B(t' ,t, p 1 ' 2 ) gives an indication of the de- 
gree of non-Markovianity of the dynamics, as expressed 
by the measure Af(A) in Eq. (17). Any variation of 
the trace distance on a finite time interval can be lower 



bounded through Eq. (10), so that one could introduce 
a lower bound to Af(A^ This is larger the more the 
effects of system-environment correlations and environ- 
mental evolution rise above the upper threshold in Fig. [I] 

As a further remark, let us notice that our analysis 
relates the (non-)Markovianity of a given evolution with 
the general correlations between the system and its envi- 
ronment, regardless of their (quantum or classical) na- 
ture. The relevant quantity B(tl ,t, p 12 ) is defined in 
terms of the total correlations Xg B (t), which include both 
classical and quantum correlations. In general, the inter- 
action does not have to build up quantum correlations 
in order to determine a reduced non-Markovian dynam- 
ics [37]. 



D(t,p l - 2 )-F(t\t,p^) 



Non-Markovian 



/ 



/ \ 




Markovian 



D(t,p'-') + F{t',t,p L 



FIG. 1. The horizontal axis represents the possible values of B 
defined in Eq. |5| at fixed times t and t' and initial reduced 
states ps(0) and p|(0), while the two marks represent the 
thresholds given by D(t, p 1 ' 2 )— F(t', t, p 1 ' 2 ) (blue region, lower 
mark) and D(t,p 1,2 ) + F(t' ,t, p 1 ' 2 ) (red region,upper mark), 
with D(t,p^ 2 ) and ,F(t', t, p 1 ' 2 ) defined in Eq. (7) and Eq. Q 
respectively. If B(t' , t, p 1 ' 2 ) lies above the upper threshold for 
some t, t',pg 2 (0), the dynamics is non-Markovian, while if it 
takes values below the lower threshold for all the t, t' , p s ' (0), 
the dynamics is Markovian. Values of B(t' ,t, p 1 ' 2 ) between 
the two thresholds are compatible with both Markovian and 
non-Markovian dynamics. 



II. EXAMPLES 

A. Transition from Markovian to non-Markovian 
dephasing 

We now consider, as our first explicit example, the 
model described in Ref. |22) . where the transition from 
Markovian to non-Markovian dephasing dynamics has 
been experimentally realized by modifying the initial 
state of the environment. The total system under in- 
vestigation consists of single photons generated by spon- 
taneous parametric down conversion and passing through 
a Fabry-Perot cavity mounted on a rotator and then 
through a quartz plate. The open-system qubit is en- 
coded in the space spanned by two orthogonal polar- 
ization states, while the frequency of the photons is 
used to encode the environment. The dephasing dynam- 
ics experienced by the polarization of light and due to 
the quartz plate can be described through the unitary 
system-environment evolution 



U(t)\X,u) = e in ^\\u) (A = H,V), 



(18) 



where |iJ),|V) denote horizontal and vertical polariza- 
tion states, n\ is the refractive index of the plate for a 
A-polarized incident photon, and \lo) stands for the envi- 
ronmental state at set frequency lo. The initial state of 
the environment is pe(0) = |'0e(O))(V , e(O)| with 



|Vb(0)) = / dw/(w)|w), 



(19) 



and the frequency distribution |/(w)| 2 is controlled via 
the tilting angle of the Fabry-Perot cavity. For any fixed 
amplitude /(w), the dynamics of polarization can be de- 
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scribed by a family of CPT maps {A(t)} t>0 such that 

A(t) ( PHH PHV ) = ( J" H k * {t)pHV ) , (20) 
\ PVH Pvv ) \ k{t)p VH p vv J 

where p\\> = (A|p s (0)|Y) (A, A' = H,V) with p 8 (0) 
the initial polarization state, and the time-dependent de- 
phasing function 



k(t) 



dw|/(w)re 



(21) 



has been introduced. Regardless of the choice of f{uS), 
the pair of initial system states that maximizes the in- 
crease of trace distance is 



1 

71 



|tf±(0)> = ~{\H)±\V)). 



(22) 



The calculation of the corresponding trace distance at 
a time t is straightforward and leads us to W(A) = 
X^(IM^')I — |fc(aj)|), where the sum is taken over the 
temporal region of extremes cij and bj > cij where \k(t)\ 
grows. The dynamics at hand is thus non-Markovian iff 



|*(t + 01 >!*(*)! 



(23) 



for some t,t' > 0. By adjusting the distribution J(lj), 
one can arrange for a transition from Markovian to non- 
Markovian open-system dynamics |22j . which can be 
characterized through the time-local master equation 



at 



Ps(t) = -ie{t) [<T z ,Ps(t)] +l{t) (cr z ps(t)a z 



-Ps(t)), 
(24) 
7(*) = 



where e(t) = — hlm.[dt lnfc(i)] and 
-iRe[9 t lnfc(t)]. 

Let us now address how our general analysis applies 
to this model. The transition from Markovian to non- 
Markovian dynamics is shown by following the evolu- 
tion of the trace distance for the pair of initial reduced 
states in Eq. (22) upon variation of the initial state of 
the environment, as set by /(w). Thus, we can focus 
on the evolution of the initial total states IV'sb(O)) — 
|V4(0)) ® \M0)) and |^(0)) = |VS(0)) ® |^ B (0)). 
The corresponding total states at time t are found using 
Eq. (18). In turn, this allows us to evaluate analytically 



all the quantities of interest introduced in Sec. |TJ which 
are defined in terms of the trace norm of operators on the 
Hilbert space of the open system. It turns out that the 
two environmental states are equal at all times, so that 
the possible increases of the trace distance can be traced 
back solely to the correlations built during the interac- 
tion between the open system and the environment. In 
particular, we find that the effects of system-environment 
correlations at time t on the subsequent evolution of the 
open system are quantified through 



B(t' ',t, p 1 ' 2 ) = \k{t + t') - k{t)k{t')\ 



(25) 



Moreover, we have that F(t' ,t, p 1 ' 2 ) = \k(t)k(t')\ and 
&D(t' ,t, p 1 ' 2 ) = \k(t + t')\-\k(i)\. Therefore, Eq. (pj) is 



indeed satisfied. These quantities can be explicitly evalu- 
ated by specifying the frequency distribution |/(w)| 2 that 
determines kit) through Eq. (21). As a first example, we 
use a Lorentzian distribution: as we shall see, this leads 
to a semigroup evolution and thus provides a natural 
benchmark. Explicitly, we take 



l/MI 2 



Slu 



7T[(W 



5uj 2 



(26) 



where w° is the central frequency and 5uj is the width of 
the Lorentzian distribution, which can be obtained using 
single photons emitted by quantum dots . Choosing a 
Lorentzian frequency distribution corresponds to taking 



k(t) 



which entails the exponential decay 
In this case, 



-8ut 



of the trace distance D(t,p ' ) = 
the coefficients of the time- local generator in Eq. (24) 
reduce to e(t) — loq/2 and j(t) = Suj/2, so that the 
dynamics of the open system is fixed by a completely 
positive semigroup [351 EH] • In addition, the exponential 
expression of k(t) means that the correlations between 
the open system and the environment have no influence 
on the evolution of the trace distance between reduced 



states, as it follows from Eq. (25). Thus, for the model at 



hand, in the semigroup regime system-environment cor- 
relations, despite being present, do not affect the dynam- 
ics of the open system at all. The exponential expression 
for k(t) leads to B(t' ,t, p 1,2 ) = 0, which in turn implies 
D(f + i,p 1 ' 2 ) = F(t',t,p^ 2 ): States p> SE (t) (j = 1,2) 
can be replaced with p J s (t) <g> psit) without modifying 
the subsequent evolution of the trace distance. 
Let us now consider the frequency distribution 



l/HI 2 - E 



,0\2 



so that k(t) 



-5(jJi)t 



O (ioj9—5ujz)t' 



/(i 



(27) 



with 



r = Ai/A2- As it can be easily checked, if the two 
Lorentzian distributions in Eq. (27) have the same cen- 



tral frequencies, = lo^i the resulting dynamics is still 
Markovian as |fe(£)| is monotonically decreasing. Inci- 
dentally, differently from the case of a single Lorentzian 
distribution, the family of CPT maps determined by 



Eq. ( 20 1 is no longer a semigroup but a divisible fam- 
ily of completely positive dynamical maps [20] . In fact, 
the coefficients of the corresponding generator in Eq. ( 24 ) 
are given by 



e(t) 



UJQ 

2 ' 



7(*) 



(28) 



which are positive at all times. Our analysis allows to 
trace the Markovianity of the dynamics back to the weak- 
ness of the system-environment correlations created by 
the interaction. In Fig. [2] (a), we plot B(t' ,t, p 1,2 ) and 
Dfap 1 ' 2 ) - Fty'^p 1 ' 2 ). Both the quantities refer to 
the two initial states of the system in Eq. ( 22 ) and are 



plotted as a function of t and t' . The effects of system- 
environment correlations on the reduced system, as quan- 
tified through B(t' ,t, p 1 ' 2 ), are always weaker than the 
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(a) (b) (c) 




FIG. 2. Plot of the effects of system-environment correlations described through the function B(t' ,t' p 1 ' 2 ) given in Eq. | |25| 
corresponding to the example discussed in Sec. |II A| for a frequency distribution as in Eq. ( |27[ ). Panel (a): we have taken 
cjJ/(5cji = ui2/8uii, Slj2/Sloi = 10, and r = 1. In both panel (b) and (c) we have taken 8ui,2 = Suj, uji/Slj = 1 and uj^/Suj — 9 
with r = 1 [t'8u> = 0.3] in panel (b) [panel (c)]. Semi-transparent surfaces represent D(t,p 1,2 ) — F(t' ,t, p 1 ' 2 ) [panel (a)] and 
D(t,p h2 ) +F(t',t,p 1 > 2 ) [panels (b) and (c)]. 



(a) (b) (c) 




FIG. 3. (a) We consider a quantum spin chain of N + 1 particles partitioned into a two-level system spin embodied by particle 
and a finite-size environment (provided by the rest of the chain). The dynamics experienced by spin can be changed 
from Markovian to strongly non-Markovian by adjusting the parameters entering the system-environment coupling model Hse 
and the inter-environment one He- Details of the form of such Hamiltonians are given in the body of the paper, (b) [(c)] 
Comparison between the function B(t' ,t, p 1 ' 2 ) (full-colorer surface) and the threshold for Markovianity [non-Markovianity] 
D(i,p 1,2 ) - F(t', t, p 1 ' 2 ) [D(t,p ia ) + F(t' ,t, p 1 ' 2 )]. (semi-transparent curves). We have used N = 8 with J /J = 1 and 
B/J=10" 2 . 



threshold embodied by D(t,p l < 2 ) - F(t' ,t, p 1 ' 2 ). As a 
consequence, the distinguishability between the states of 
the open system cannot increase and a Markovian dy- 
namics is induced. 

A different situation occurs if the two Lorentzian dis- 



tributions in Eq. (27) have the same width but dif- 
ferent central frequencies. In this case, the dynamics 
is non-Markovian, i.e. \k(t)\ is a non-monotonic func- 
tion of time. As shown in Fig. [2] (b), the effects of 
system-environment correlations are now stronger than 
the upper threshold D(i,p 1 ' 2 ) + F(f ,t, p 1 ' 2 ). Basically, 
there are times t such that the two reduced states, Pg(t) 
and pi(t) are very similar, so that the upper threshold 
D(t, p 1 ' 2 ) + F(t' , t, p 1,2 ) is small, while the correlations of 
the two total states, XsbW an d are stm differ- 

ent. Consequently, the trace distance at subsequent times 



t + t' increase. In this case, system-environment corre- 
lations do induce a non-Markovian dynamics. Finally, 
by taking r in Eq. (271 from zero to any non-zero value, 



we have a transition from a Markovian dynamics, more 
precisely a semigroup dynamics, to a non-Markovian dy- 
namics. In Fig. [2] (c) one can see how this is reflected in 
the behavior of B(t' , t, p 1,2 )- In fact, the latter turns from 
being identically zero for r — to increasing above the 
upper threshold D(t, p 1,2 ) + F(t', t, p 1 ' 2 ), thus implying a 
non-Markovian dynamics. 



B. Non-Markovianity in a spin-chain system 

As a second instance, we consider the case of a single 
qubit attached to a quantum spin chain of N spin- 1/2 



7 



particles, along the lines of the studies reported in Ref. [5J 
[34] and as illustrated in Fig. 3] (a). The open system is 
embodied by spin 0, while the environment consists of 
particles 1 — > N. The overall N + 1 system are mutually 
coupled via an XX model and subjected to a transverse 
magnetic field. Assuming units such that h = 1, the 
corresponding Hamiltonian model is H—Hs e + He with 



H. 



N-l 



N 



H E = -2 J J2 (KK + i + ZZK+l) " 25 E *n, 



(29) 



n=l 



n=l 



where cr\ is the /c-Pauli matrix (k=x, y, z) for particle n, 
B is the amplitude of the magnetic field affecting S and J 
(Jo) is the inter-environment (system-environment) cou- 
pling strength. The underlying assumption is that the 
free evolutions of S and E are identical, thus allowing the 
passage to the interaction picture without the introduc- 
tion of time-dependent coefficients. The non-Markovian 
evolution experienced by spin was characterized fully 
in Ref. [5], where it was found that, for interaction times 
that are within the recurrence time of the system (when 
any information propagating across the chain returns to 
the open system after reaching the end of the chain), 
there is a working point defined by (Jo/J, B / J) at which 
the measure of non-Markovianity J\f is null. As the op- 
timization inherent in the definition of such measure is 
achieved for system states lying on the equatorial plane 
of the Bloch sphere [S]) we consider the input states 
/°s(0) = IV's (0)}(V , s (0)|, while the environment is ini- 
tialized in Pi >2 (0) = pfni = «'ili 1 |0)i(0|. In order to pro- 
vide a physically significant example that is nevertheless 
able to show clearly the features that we are interested 
in, we solved fully the problem embodied by an environ- 
ment of TV = 8 spins with J /J — 1 and B/J = 10~ 2 . 
At variance with the previous example, the system lacks 
of an analytically amenable solution (the expressions for 
the trace distance and the quantities introduced in the 
Section [I] are too involved to be reported here) but al- 
lows for a handy numerical analysis. The results are 
shown in Fig. 3] (b) and (c), where B(t,t',p 1 ' 2 ) is com- 
pared to the thresholds D(t, p 1 ' 2 ) ± F(t, t', p 1 ' 2 ) within a 
broad range of values for t and if. Clearly, besides the 
existence of a range of values of (t, if) where the dynam- 
ics is expectedly Markovian, B(t, t' , p 12 ) soon trespasses 
the thresholds for non-Markovianity [cf. Fig. [3](c)]. In 
particular, this is the case for values of (t, if) such that 
J • (t + if) < 3, which guarantee that the correspond- 
ing evolutions occur well within the recurrence times of 
the spin chain and any non-Markovian effect is due to the 
intrinsic features of the interaction rather than the finitc- 
ness of the environment. Interestingly enough, the gap 
between the lower and upper threshold identified above 
may disappear, in this example. Values of t exist at which 
D(i,p 1:2 ) ~ F(t,t' , p 1 ' 2 ), thus making the gap between 
upper and lower threshold effectively null. 



III. CONCLUSIONS 



We have shown how the dynamical correlations be- 
tween an open quantum system and its environment in- 
fluence the nature of the open system's dynamics, as 
far as non-Markovianity is concerned. Our analysis re- 
lies on the definition of non-Markovianity given in terms 
of the evolution of the trace distance between reduced 
states [T§] . 

We have introduced the quantity F{if ,t, p 1,2 ) [cf. 
Eq. Q] that describes the evolution of the trace dis- 
tance under the CPT maps that would connect states at 
different times if correlations and environmental changes 
could be neglected at any time. In a complementary way, 
we have defined the quantity B(t! ,t, p 12 ) [see Eq. ([5|] 
that measures, by means of the trace norm, the ef- 
fects of system-environment correlations and environ- 
mental evolution due to the interaction up to a time t 
on the subsequent dynamics of the open system. These 
two quantities allow to introduce an upper and a lower 
bound to the variation of the trace distance between re- 
duced states on finite time intervals, as quantified by 
Eq. ( 10 ). We have thus been able to conclude that if the 



effects of correlations and environmental evolution are 
below a first threshold the resulting reduced dynamics 
is certainly Markovian. Despite being necessary, system- 
environment correlations and changes in the environment 
are not a priori sufficient to induce an increase of the 
trace distance. On the other hand, if their effects exceed 
a second threshold, a non-Markovian reduced dynamics 
is surely induced. 

The general analysis has been applied to the model ex- 
ploited in |22j to experimentally detect the transition be- 
tween Markovian and non-Markovian dynamics. We have 
shown how such transition can be explained in terms of 
the different effects of system-environment correlations. 
By properly varying the initial state of the environment, 
one can describe both a semigroup dynamics, in which 
system-environment correlations do not affect at all the 
evolution of the reduced state's distinguishability, and 
non-Markovian dynamics, in which system-environment 
correlations strongly influence the dynamics of the open 
system. 

Our results will be useful also with respect to different 
approaches to non-Markovianity relying on other proper- 
ties of the dynamical maps, since they show in full gen- 
erality to what extent system-environment correlations 
can be compatible with a contraction of the trace dis- 
tance. In particular, this could help to further under- 
stand the connection between correlations in the total 
state and breaking of divisibility of the completely pos- 
itive dynamical maps [SJ [T^]. In addition, our results 
could provide further insights into microscopic deriva- 
tions of reduced dynamics, in order to clarify the role 
of system-environment correlations and changes in the 
environmental state 140 . 
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